Influence of long-range interactions on the critical behavior of the Ising model 
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We study the ferromagnetic Ising model with long-range interactions in two dimensions. We first 
present results of a Monte Carlo study which shows that the long-range interactions dominate over 
the short-range ones in the intermediate regime of interaction range. Based on a renormalization 
group analysis, we propose a way of computing the influence of the long-range interactions as a 
dimensional change. 
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Introduction — In recent years there has been a lot of 
interest in the statistical physics of classical and quantum 
systems with long-range interactions, for a review see [l[. 
The role of quasistationary states and ergodicity break- 
ing in long-range interacting systems was investigated in 
1 and i. In [J] the approaching to equilibrium for long- 
range quantum systems was examined and there has been 
a lot of enthusiasm in investigating the entanglement en- 
tropy in long-range spin chains [5|-|7[. Very recently an 
experiment was conducted on a quantum system with 
tunable long-range interactions Q. 

In the present study wc focus on the Ising model which 
is probably the most studied model in statistical me- 
chanics, especially in the context of critical phenomena. 
Most of the studies about the Ising model are concen- 
trated around the short-range case which is exactly solv- 
able in one and two dimensions Q. In three dimen- 
sions the problem was perturbatively studied using the 
e-expansion technique [1 01 ] of the renormalization group 
(RG) combined with the Borel resummation of the per- 
turbation series, see [TTJ] and references therein. Most 
recently the problem was revisited by using conformal 
bootstrap technique [13]. Although now there are lit- 
tle unknown facts around short-range Ising model the 
long-range Ising model is still the subject of many con- 
tradicting theoretical and numerical studies. We define 
the long-range Ising model as 
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where the sum is over all pairs of spins of a d dimensional 
system and J > . In [13[ the e-expansion technique was 
applied to the above problem shortly after the introduc- 
tion of the method. Three regimes were discovered: (a) 
the classical regime < a < d/2 with mean field criti- 
cal exponents; (b) the intermediate regime d/2 < a < 2, 
where the exponents are functions of a and (c) the short- 
range regime a > 2 where the exponents are the same as 
in the short-range Ising model. The conjectures around 



the first regime are already proved in [l4j and the results 
of the third regime are widely accepted. The interme- 
diate regime has been the subject of many controver- 
sies in the last forty years. Based on [l3| the r\ expo- 
nent in this regime follows the classical regime's formula; 
i.e. i] = 2 — a and apparently it has a jump at a = 2 
to the short-range value rj SR > in two and three di- 
mensions. Shortly after, Sak [lj| argued that there is no 
such discontinuity in the exponent because if one looks at 
the long-range interaction as a perturbation of the short- 
range Ising, cr = 2, one can discover that the short-range 
Ising exponents should be extended up to a = 2 — r/ SR and 
so there is no discontinuity in the critical exponents. Al- 
though some other forms of RG also appeared [l6j] in the 
last forty years, Sak's argument has been widely accepted 
[13] ■ An especially interesting numerical simulation done 
by Blote and Luijten [3| completely ruled out any jump 
in the value of rj. The motivation of our work comes from 
a recent numerical work done by one of us [I3| in which 
we improved the algorithm, analyzed much bigger sizes 
and concluded that neither Fisher et al. procedure nor 
Sak's machinery fit with the numerics, in particular in 
the intermediate regime and close to the boundary with 
the short-range regime. 

In the present study, we will concentrate on two aspects 
of this problem. First, we will compare the long-range be- 
havior with the short-range one in two dimensions. Wc 
provide numerical evidences that the long-range interac- 
tions dominate for a < 2. Next, we will propose another 
way to compute the r/ exponent. The main idea is to 
make a correspondence between 

M = Jd d x Q|V*/ 2 Sb (x)| 2 + A/ Mx) |^ ; (2) 
with e' = 2o~ — d and 
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with e = 4 — D where we introduced fi to make A di- 
mensionless. The first expression A\ is a formal way 
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of writing in real space a model with long-range inter- 
actions. The second expression A2 is the usual way of 
expressing a <j) A model for a D dimensional theory. For 
(7 — 2, and with the condition e = e', the e-expansion of 
both models, A\ and A2 will be the same apart from a 
term proportional to 5a = 2 — a. Thus the computation 
will be done from the model A2 with D = 4 + d — 2a. 
The deviation of a from 2 is replaced by a deformation 
of the dimension from d to D. 

Long-range versus short-range — In it was ob- 

served that the behavior of the model with long-range 
interactions and for a < 2 is different from what is ex- 
pected for the short-range model. Then we must worry 
about the relevance of the long-range compared to the 
short-range. If we start from a short-range model and 
consider the addition of long-range term g SiSj/rf^ a 
as a small perturbation, then a simple dimensional argu- 
ment predicts the relevance_afthe perturbation as a func- 
tion of the dimension of g 
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17| . Since for large distances 
with 77 SR = j in two dimen- 
sions, we expect that the dimension of g is 2 — r\ 

SR . 

Then the result is that the perturbation is relevant for 
a < 2 — 77 S R and irrelevant otherwise. We will now test 
this argument. We consider the case of the perturba- 
tion as a function of a and g. We will use the magnetic 
cumulant defined as: 
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where K = j3J. For each value of g, we consider the 
quantity B c (g, L, L') which corresponds to the crossing 
of B(g,L,K) and B(g, L' , K) as a function of K, By 
choosing a set of increasing values L and L' not too far 
apart, we can determine for each pair the value of K 
which corresponds to the crossing and B c (g, L, L') is ex- 
pected to converge towards a finite limit for L — > +00. In 
the following, we will always consider L' = 2L and then 
we will just denote the crossing value by B(g, L). In fl9| . 
it was determined that the corresponding quantity for the 
long-range interaction model, which can be considered as 
the limit g — > 00, converges to a value smaller than the 
one of the short-range model for a < 2. In Fig. [TJ we 
present the measured values of B(g,L) for a = 1.6 and 
a = 1.8. For the first case, we observe a clear tendency 
for B(g, L) to converge towards the same limit as the LR 
model (which is shown as a dotted line) and this for all 
the values of g in the range g = 0.01uptog=l. For 
the second case, the situation is less clear. For the small 
values of g < 0.1, it seems first that B(g,L) converges to- 
wards the model with short-range interactions. For larger 
values of the perturbation, we just observe that B{g 1 L) 
increases with the size. While it can be assumed that this 
just corresponds to the flow towards the value for the SR 
model, one can also invoke the effect of strong finite size 
corrections. 

In fact, since we are considering a case in which there 
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FIG. 1: (Color online) B(g,L) vs. L for a — 1.6 (top panel) 
and a — 1.8 (bottom panel). The SR results are depicted in 
red and have the largest B(g, L). The dotted line corresponds 
to the extrapolated value for the LR model (see JT?" 



is both a flow towards either the LR model or the SR 
model and very strong finite size effects, it is difficult to 
know which one is the dominant effect. We then adopt 
another strategy. We will look in the following to the 
quantity defined as 



X(g,L) 



B(g, L) — B hR (L) 
B SR (L)-B LR (L) 
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This quantity is defined such that it takes a value be- 
tween and 1. If B(g,L) flows towards the SR point, 
then X(g,L) goes to 1. On the contrary, if it flows to- 
wards the LR point, then it goes to 0. We then expect the 
crossover to be controlled by a crossover parameter g/\t\^ 
with t the reduced temperature and (f> = i^A CT which de- 
fines A CT . On a finite lattice of linear size L and at the 
critical point, this becomes gL Aa with the correspon- 
dence t~ v ~ £ ~ L. According to the naive dimensional 
analysis made above, A CT = 2 — a — r/ SR . In Fig. [21 we 
show a plot of X(g,L) vs. the crossover parameter for 
various values of g in [0.01, 1] and a. For each value of 
tr, we determined a single parameter A a which allows to 
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FIG. 2: (Color online) X(g,L) vs. gL A<r . The straight dashed 
line corresponds to (gL A °)~ ' 75 , while the other dashed curve 
is for exp(— gL A "). 



part, with respect to the renormalized coupling and held 

as 
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where e' = 2a — d, s = Z —1 ' 2 Sb and Ao = A/i £ Z\ with 
Sb and Ao as bare parameters. The scale p is introduced 
because we want to do the expansion around the massless 
theory [21| . The renormalization conditions are 
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make a scaling for all the values of g on a single curve. 
The values that we obtained are reported in the cap- 
tion of the figure. It is quite remarkable that the curves 
for all values of a collapse on a single curve. We ob- 
tain that for gL Aa <^_ 1, the curves follow an exponential, 
i.e. X(g,L) ~ exp(-gi A<T ). For pL A<T > 1, the curves 
behave as a power law i.e. X(g 7 L) ~ (gL A <*)~ a with 
a ~ 0.75. 

A second fact is that the value of A a does not follow the 
prediction obtained from the naive dimensional analysis. 
While for small values of a, the correspondence between 
the measured crossover exponent and the predicted one 
is acceptable, this is clearly not the case for larger values 
of a. And in particular, this exponent does not cancel 
at a = 2 — 7] SR . Note also that the precision on this 
exponent is not very good for larger values of a since in 
that case the denominator of X(g, L) becomes small and 
will cancel in the large size limit for a = 2. 

The conclusion of this analysis is that we observe a 
clear signal for a crossover between the short-range inter- 
action model and the long-range interaction model. This 
crossover seems to be present in all the range that we 
can consider 1.2 < a < 1.9. Of course, such a crossover 
is expected for small values of cr, i.e. a < 2 — ?y SR = 1.75. 
We observed that in fact this crossover between the SR 
model towards the LR model remains present even for 
larger values of cr, presumably up to a = 2.0. 

Renormalization group approach — In this section we 
propose a new way of doing RG analysis for long-range 
Ising model. Although our analysis shares some similari- 
ties with the work of Yamazaki the final results are more 
general [l6|. We implement our RG analysis around the 
critical point [2ll . 1 2 21 ] , then we can avoid calculating more 
complicated integrals. Based on the arguments of the last 
section one can write the Lagrangian of the long-range 
Ising model, forgetting about the irrelevant short-range 
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of the Fig. [3] one can write 
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The integrals are infrared divergent for e' > and need 
to be calculated by analytical continuation from the con- 
vergent region. This is the same situation as in the usual 
e-expansion in the short-range Ising model. Although the 
integrals are complicated, they can be calculated using 
the formulas in [21[, and after using the renormalization 
conditions we will have 
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A 2 



-] r[ 



Z\ = 1 



3(r ( 47r )d r[Mr3£] r 3 [f] 

3A r 2 [^] r[e'/2] 
2(47r) rf / 2 r[d-cr] r 2 [a/2]' 



(12) 
(13) 



The very important point is that if we expand Z with 
respect to e' there will not be any pole and one cannot get 
sensible contribution to the critical exponent of the field 
s(x). However, since the integrals are infrared divergent, 
the right way to get a sensible perturbation theory is to 
expand Z first around a = 2 and then around e' = 0. 
Since we have two parameters dimension d and cr; and 
they can be changed independently, one can first consider 
having 5a = 2 — a small and then do the perturbation 
theory with respect to e'. After expanding Z and Z g with 
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respect to 5a and then e' we get 



A 



1 



12(47r) d e' 
3A 

+ ( 47r )d/2 1 



1 

2~ 



127-13' 



(5(7- 



0(A 3 ), 
(14) 



7 [1 + (1 -7)^ + ••■]+ 0(A 2 ). (15) 



with 7 = 0.5772... as the Euler-Mascheroni constant. Us- 
ing the Callan-Symanzik equation which states that the 
derivative of the bare quantities with respect to \x is zero, 
one can get the beta functions as 
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To derive the above formula we first use the equation (15) 
to get a relation between A, Ao and fi. Using the above 
beta functions at the critical point where /3(A*) = one 
can easily get the correction to the mean field value of 
the critical exponent Srj = r/ — (2 — a) as 
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Based on our prescription it is obvious that in the e'- 
cxpansion of the 77 exponent the zeroth order terms of 
5a expansion will be the same as the e-expansion of the 
short-range Ising model but with e' = 2a — d instead of 
e = 4 — D. So in principle close to the a = 2 we will have 

r 1 = 2-a+^-e' 2 + --- + 0(5a), (19) 
54 

where dots represent the higher order terms of the e'- 
cxpansion. Since in our proposal of doing the RG we 
first expand all the contributions around a = 2 and then 
around e' = we expect that the dots in the formula (TT9")) 
are exactly the same as in the short-range Ising model 
with e' instead of e. The above expansion suggest that for 
5a small one can argue that the critical exponent of the 
long-range Ising model in d dimensions is approximately 
the same as the critical exponent of D = 4 + d — 2a short- 
range Ising model. For the short-range Ising model 5rj is 
known up to e 5 for various dimensions [23, 24 1. The first 
correction to this value comes from the second and third 
terms of the equation ([T5)l which are both negative. If 
the higher order terms, (5a) n with n > 2, do not change 
the sign of the contribution to 77 one can conclude that 
the critical exponent 77 of the short-range Ising model 
in D = 4 + d — 2a gives an upper bound for the 8rj of 
the long-range Ising model in d dimension. Of course this 
conjecture needs to be checked by calculating higher loop 
corrections to the critical exponents. Based on the above 
arguments we compared in Fig. U the r\ coming from the 
numerical calculations for the long-range Ising model in 



two dimensions with the results coming from the five loop 
calculation of D = 6 — 2a dimensional short-range Ising 
model. The results are well comparable in the region 
1.75 < a < 2 and as we argued for the smaller values of 
a the actual values lie below our approximation. 





FIG. 3: The relevant Feynman diagrams in RG calculation of 
the beta function and wave function normalization. 




FIG. 4: Data depicted in red are the result of the simulation in 
[Tgj ] . Green and blue data are coming from our approximation 
using the results of and [24|] respectively. The dotted line 
is Sak's result El. 



Conclusions — In this letter we provided further nu- 
merical evidences that the long-range interaction have an 
influence on the critical behavior of the Ising model for 



a < 2, in contrast with previous RG studies [13j, [TpJ. We 
proposed a way to compute the influence on 77 of a devi- 
ation from a = 2 based on renormalization group ideas. 
The main idea is the double expansion with respect to 
5a = 2— a and then e' = 2<t— d in a way that we get a non- 
trivial contribution to the wave function renormalization. 
Our analysis shows that close to a = 2 one can approx- 
imate the 77 exponent of the d dimensional long-range 
Ising model with the same exponent of d = 4 + D — 2a 
dimensional short-range Ising model. Our results are in 
excellent agreement with the numerical results (l9j . 
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